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ABSTRACT 

The Northern Sky VariabiUty Survey (NSVS) is a temporal record of the sky over the optical mag- 
nitude range from 8 to 15.5. It was conducted in the course of the first generation Robotic Optical 
Transient Search Experiment (ROTSE-I) using a robotic system of four co-mounted unfiltered telephoto 
lenses equipped with CCD cameras. The survey was conducted from Los Alamos, NM, and primarily 
covers the entire northern sky. Some data in southern fields between declinations 0° and —38° is also 
available, although with fewer epochs and noticeably lesser quality. The NSVS contains light curves for 
approximately 14 million objects. With a one year baseline and typically 100-500 measurements per 
object, the NSVS is the most extensive record of stellar variability across the bright sky available today. 
In a median field, bright unsaturated stars attain a point to point photometric scatter of '^0.02 mag and 
position errors within 2". At Galactic latitudes |5| < 20° the data quality is limited by severe blending 
due to ~14" pixel size. We present basic characteristics of the data set and describe data collection, 
analysis, and distribution. All NSVS photometric measurements are available for on-line public access 
from the Sky Database for Objects in Time-Domain (SkyDOT; http://skydot.lanl.gov) at LANL. Copies 
of the full survey photometry may also be requested on tape. 

Subject headings: catalogs - surveys - astronomical data bases - stars: general, variables 



L INTRODUCTION 

An amazing fact about modern astronomy is that the 
global time variability of the optical sky is largely un- 
explored for objects fainter than those observable with 
the naked eye (Paczyhski 1997). As a result, the exist- 
ing samples of known variables are quite incomplete. The 
commonly accepted standard catalog of variable stars is 
the General Catalog of Variable Stars (GCVS; Kholopov 
1998). However, the GCVS catalog was compiled from 
a multitude of heterogeneous observational sources — in 
many cases based on analysis of photographic plates — 
and does not present any light curves. Spatial coverage 
in the GCVS is strikingly patchy at magnitude 10 and 
fainter. With the data currently available for most stars, 
we cannot answer the simple question of whether or not 
the star is variable. This is unfortunate because tempo- 
ral fiux changes not only carry useful physical information 
about the star but they are of concern for experiments 
where variability could degrade the quality of the compar- 
ison star grid (e.g. Space Interferometry Mission, Frink et 
al. 2001). 

The reasonably sensitive, inexpensive CCDs that have 
become available in the last few years, coupled with 
the improved affordable data processing capabilities, have 
opened new windows for discovery in astrophysics. The 
large, uniform samples being collected with CCD imagers 
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are enabling the temporal study of objects in new detail. 
Microlensing surveys, for example, have shown that mas- 
sive photometric monitoring programs can return numer- 
ous scientific results, often unrelated to the original goal 
(Paczyhski 2000a, Ferlet, MaiUard & Raban, eds. 1997). 
The added value of large number statistics and good sky 
coverage is also evident in catalogs from digitized POSS 
plates (Djorgovski et al. 2001) and multicolor surveys like 
SDSS (Stoughton et al. 2002, Abazjian et al. 2003) and 
2MASS (Skrutskie et al. 1997) when apphed to rare ob- 
jects such as high-z quasars and galaxies or brown dwarfs. 

A new generation of small robotic sky patrol instru- 
ments is making all sky temporal monitoring of point 
sources possible (Paczyhski 2000b, Chen, Lemme & 
Paczyhski, eds. 2001). New results from these variabil- 
ity studies are contributing to our understanding of stellar 
evolution and Galactic structure. Resulting improvements 
in the local distance scale and discoveries of supernovae en- 
able more accurate estimates of cosmological parameters 
and stellar ages. A substantial contribution to the under- 
standing of gamma-ray bursts (GRBs) has already been 
made by robotic follow-up telescopes of modest size like 
ROTSE (Akerlof et al. 2000b, Kehoe et al. 2001) and 
LOTIS (Park et al. 2002). Development of autonomous 
systems searching for optical flashes in real time will en- 
able monitoring of a variety of fast and rare phenomena 
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including the onset of optical emission from GRBs. The 
RAPTOR system (Vestrancl et al. 2002) is a stereoscopic 
sky monitoring system that is designed to find such flashes 
and provide instant notification and response while events 
arc occurring. In the context of extra-solar planets, a 
shallow but very large area time domain survey with high 
cadence and sub-percent photometry currently offers the 
best prospect of discovering bright systems with transiting 
Jovian planets that are suitable for detailed studies using 
high resolution spectroscopy (Horne 2002, Charbonneau 
2002). 

To date, the largest volumes of data on time variabil- 
ity are those collected by microlensing searches, contain- 
ing typically 10^ light curves, each with a few hundred 
individual photometric measurements spread over several 
years (Paczyhski 2001, Brunner et al. 2001). Microlensing 
surveys can be considered medium angle surveys limited 
to specific areas of interest, namely the Galactic Bulge and 
several galaxies of the Local Group. In the category of very 
wide field surveys, the All Sky Automated Survey (ASAS, 
Pojmahski 1997) has made a substantial fraction of the 
data available in the public domain and returned ~7,000 
variable stars brighter than 15 mag discovered primarily 
in the 6'' quadrant and additional scattered fields of the 
southern hemisphere (Pojmahski 2000, 2002). The search 
for transits by extra-solar Jupiters resulted in recent pro- 
liferation of large area photometric monitoring projects 
(Horne 2003), however little data has been published to 
date. 

This paper marks the first release of data from the 
Northern Sky Variability Survey (NSVS), a CCD-based 
synoptic survey covering the entire sky north of the decli- 
nation S = —38°. This data release provides light curves 
for ~14 million objects down to V ~15.5 magnitude, with 
hundreds of repeated observations spanning one full year. 
NSVS photometry is a great improvement over photo- 
graphic work despite disadvantages of using a single un- 
filtered photometric band. Most importantly, all obser- 
vations were collected with the same instrument resulting 
in uniform data quality limited largely by crowding near 
the Galactic plane. We are working toward making the 
NSVS the best possible tool for studies of stellar variabil- 
ity and Galactic structure using bright stars. Akerlof et 
al. (2000a) published a preliminary variability analysis 
covering 9 out of 161 survey tiles and 3 months of observ- 
ing time. While we repeat parts of their discussion of the 
observing system and data reduction, there are significant 
differences in data processing between this work and Ak- 
erlof et al. (2000a). 

The outline of the paper is as follows. Section 2 de- 
scribes the telescope, cameras and the process of data col- 
lection. Section 3 presents the details of data reductions 
and photometry followed by the discussion of survey qual- 
ity and coverage. In Section 4 we describe public access to 
the data and some technicalities of the data products. Sec- 
tion 5 concludes with the summary and future prospects 
for NSVS. 

2. INSTRUMENTS AND DATA ACQUISITION 
Numerical data for this section is summarized in Table 1. 

2.1. ROTSE-I Robotic Telescope 



All data in the present NSVS datasct were collected by 
the first generation Robotic Optical Transient Search Ex- 
periment (ROTSE-I). The primary goal of that experiment 
was prompt response to GRB triggers from satellites in or- 
der to measure the early light curves of GRB optical coun- 
terparts (Akerlof et al. 2000b, Akerlof et al. 1999, Kehoe 
et al. 2001). The normal operation of the ROTSE-I instru- 
ment was completely automatic requiring only periodic 
maintenance. The telescope consisted of four Canon 200 
mm lenses with f/1.8 focal ratio, each covering 8.2° x 8.2° 
for a total field of view about 16° x 16°. All four optical el- 
ements were carried by a single rapidly slewing mount and 
designated with the symbols a through d. The telescope 
and mount were located on the roof of a military surplus 
enclosure, which housed the instrument control comput- 
ers, and protected by a clamshell during the day or in bad 
weather. The instrument was located at the Los Alamos 
National Laboratory (W106°15'13", N35°52'9"), just out- 
side the town of Los Alamos, New Mexico. Light pollution 
at this site, although detectable, had only a minor impact 
on final survey photometry. Given the very dry local cli- 
mate and 85% fraction of useful nights it is well suited for 
survey astronomy. Significant cloud cover is confined to 
the monsoon season lasting from July to mid September. 

2.2. Imaging Cameras 

Each of the four Canon lenses was equipped with a 
thermo-electrically (TE) cooled AP-10 camera, which em- 
ploys a Thomson TH7899M CCD. The 2kx2k chip format 
covers an 8.2° x 8.2° degree field of view with 14.4" pix- 
els. The spatial resolution of the system was limited by 
instrumental seeing. The Canon lenses delivered a typical 
point spread function (PSF) with full width at half max- 
imum (FWHM) of ~20" and therefore marginally under- 
sampled images of point sources. The ROTSE-I telescope 
was operated without any filters so the spectral response 
is primarily limited by the sensitivity of the CCD, result- 
ing in a very broad optical band from 450 nm to 1000 
nm which covers the photometric bands from mid B to 
mid I. The quantum efficiency of the front illuminated 
thick CCD chips makes the effective band most compara- 
ble to the Johnson R band. To optimize the readout speed 
for GRB response measurements, the images were read in 
14-bit mode. There is no loss of information however, be- 
cause of the relatively narrow dynamic range of the AP-10 
cameras. The typical gain setting was ^8 e^ADU^^. Im- 
ages are sky background limited primarily because of large 
pixel size. The limiting V magnitude of the faintest stars 
recorded in 80-second exposures was typically 14.5-15.5. 
Saturation occurred at 10-10.5 mag in normal exposures 
(80 seconds) and at about 8 mag in bright time expo- 
sures (20 seconds). Vignetting in the lenses is very signif- 
icant and amounts to about 40% loss of sensitivity near 
the corners of the CCDs. This effect is very stable and 
easily corrected by flat fielding procedure. The shutters 
in the AP-10 cameras did not perform according to their 
specification. Especially in cold conditions, shutters did 
not operate smoothly and caused "anomalous vignetting" 
near the frame edges in some of the images. Photomet- 
ric corrections explained in Section 3.2 remove this and 
other effects. Dark frames and small scale flat field fea- 
tures were stable over a few days to a week. A few bad 
columns in the CCDs did not affect the overall quality of 
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the data set. Cameras a-d did not perform equally. The 
slightly lower photometric quality of camera d can be seen 
in survey statistics presented in Section 4.1. The loss of 
observing time due to temporary failure of camera c is also 
visible. 

2.3. Observing Protocol 

Despite the fact that the primary goal of the ROTSE-I 
project was rapid response to GRB triggers and not sky 
patrols, almost all observing time was actually spent in 
the latter mode. An accessible GRB position would be 
posted by the GCN network approximately once every 10 
days. Upon receipt of the coordinates, the ROTSE-I sys- 
tem would abort the current patrol activity and imme- 
diately start observing the field around the position for 
approximately one hour of imaging. 

At the beginning of each night, about 12 dark frames 
were collected for calibration purposes. No special flat 
field exposures were made (Section 3.1.1). The large com- 
bined field of view delivered by the ROTSE-I system re- 
quires only 206 tiles to cover the entire sky, with 161 tiles 
observable from Los Alamos. The list of fields with ele- 
vation above 20° was prepared by startup scripts. During 
normal execution of a patrol, for each of those fields two 
80 second exposures would be taken, separated only by 
the 1.5 minute duty cycle. Reduced exposure time of 20 
seconds was used in bright Moon light conditions (^30% 
of data) . No frames were taken when the moon was closer 
than 12° from the field center. Time keeping accurate to 
20 ms was implemented using Network Time Protocol. On 
a good night, it was possible to cover the entire local sky 
(~ 10,000 deg^) twice. Paired observations are useful for 
detecting variability and aperiodic transients. They also 
provide a handle on spurious detections due to man-made 
space objects, cosmic rays, hot pixels and other effects. 

The position angle of the cameras was fixed at PA=0° 
in all fields except for near-polar region where the con- 
trol software allowed PA=180°. In the latter case, fields 
normally assigned to cameras a and b would be imaged 
by cameras c and d respectively. Such observations are 
flagged appropriately and excluded from parts of the anal- 
ysis (Section 3.2.3). 

3. DATA PROCESSING 

Analysis of the data presented in this release was con- 
ducted off-line on archival ROTSE-I images. Here we 
briefly discuss the data reduction pipeline, schematically 
shown in Figure 1, which was employed to analyze that 
data. 

3.1. Image Reductions 

3.1.1. Basic Frame Corrections 

Between August 1997 and December 2001, ROTSE-I 
collected 7 TB of image data, however the performance 
of the system was not optimal in the first few months of 
the project and near the end of its lifetime. To build the 
NSVS, we selected observations covering 1 full year be- 
tween April 1999 and March 2000, when the system de- 
livered the best overall data quality. This limited the raw 
data set to 225,000 images (~2 TB). The system would 
automatically prepare a median dark frame for each expo- 
sure time using all dark images collected on a particular 



observing night. Dark subtraction removes a small frac- 
tion of pixels (■d%) with high dark current rates. Flat 
field frames were obtained from a median of all individual 
patrol images made during a given night. This is possi- 
ble with a large number of independent fields (^80) and 
statistics limited by sky noise due to large pixel size. Stel- 
lar profiles are completely removed by the procedure. We 
found that shutter problems (Section 2.2) affected some of 
the flat field images. Therefore, we visually evaluated all 
flat fleld frames and their ratios with frames made on a few 
other nights to select calibration sets with consistent large 
scale properties. As a result, it was possible to correct the 
vast majority of the frames using good calibration frames 
from the same night or the night before. 



3.1.2. Source Extraction 

The corrected images are passed to SExtractor software 
(Berlin & Arnouts 1996) which reduces them to object 
lists. The choice of this source extraction package was mo- 
tivated by the undersampling of stellar profiles, the signifi- 
cant gradients of the PSF shape in very wide field images, 
and overall reduction speed. SExtractor was optimized 
for reduction of images in galaxy surveys, but it is known 
to perform well in moderately crowded stellar fields. In 
order to optimize sensitivity, our images are filtered be- 
fore object detection with a Gaussian kernel employing a 
FWHM of 2.5 pixels and requiring a minimum of 5 con- 
nected pixels in an object. These basic detections are fur- 
ther thresholded by the software and an attempt is made 
to break up blended objects. We use SExtractor aperture 
magnitudes calculated with the 5 pixel (72") aperture di- 
ameter. Since the ROTSE-I images are almost completely 
dominated by stars, we store only a small fraction of in- 
formation available for each object: position, magnitude, 
magnitude error and processing flags. The observed errors 
at the bright end of the magnitude range are larger than 
predicted by simple photon noise. Such discrepancies are 
common for CGD measurements and are typically gener- 
ated by residual systematic effects of flat field errors, thin 
clouds, PSF variations and sampling. In order to account 
for these systematic errors, we had to add, in quadrature, 
a 0.01 mag contribution to the formal error bars. 

3.1.3. Blending 

SExtractor does not perform PSF photometry, and in 
general it is unable to deblend light distributions with- 
out a saddle point. This sets the distance limit of about 
3.0 pixels for separation of stellar blends. We found that 
the default parameters of the deblender were very conser- 
vative, resulting in very large patches of the sky at low 
Galactic latitudes being assigned to the same objects. Af- 
ter some experimentation we were able to partially con- 
trol this process, however, there are still cases when tight 
groups of several objects with merging wings are consid- 
ered to be a single object. Typically, such aggregates ex- 
tend up to ten pixels across with a bright object in the 
middle, but occasionally in dense parts of the Milky Way, 
they can be up to 30 pixels wide. This makes the com- 
pleteness of the NSVS at |6| < 20° depend strongly on 
stellar number density at small spatial scales of around 
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~ 3'. The exact assignment of pixels to objects may dif- 
fer from frame to frame and therefore some measiirements 
lost due to severe blending may still be present among "or- 
phaned" measurements, unidentified with any of the light 
curves (Section 3.2.3). The term "orphan" used here in the 
context of a photometric detection should not be confused 
with orphan GRBs and other optical transients. 

3.1.4. Astrometric and Photometric Matching 

Initial calibration of object lists consists of transforma- 
tion of instrumental positions to celestial coordinates and 
conversion of raw instrumental magnitudes to a photomet- 
ric system that can be understood by users. In wide field 
imaging this can be done on an image by image basis. 
Each ROTSE-I field covers 64 deg^ and contains on aver- 
age about 1500 stars from the Tycho catalog (Hog 1998). 
The Tycho catalog, a product of the Hipparcos/ Tycho mis- 
sion, provides very accurate astrometry and two-color {B 
and V) photometry. Most Tycho stars are fainter than the 
10 mag saturation limit of the ROTSE-I observations. The 
astrometric matching is performed in the detector plane af- 
ter deprojecting the corresponding part of the catalog us- 
ing a canonical gnomonic projection (Calabretta & Greisen 
2002). Using an approximate mount position for a given 
image, the first set of roughly 30 Tycho stars can be identi- 
fied with the triangle algorithm. This first order transfor- 
mation is used to further match a few hundred bright, but 
unsaturated, Tycho stars. A third order polynomial warp 
adequately describes the transformation between the ob- 
served and catalog positions of stars in the detector plane. 
Finally, transformed (x, y) positions are converted back 
to (q;2ooo; <520oo) using the inverse of the initial gnomonic 
projection. 

Photometric calibration is somewhat complicated by the 
very wide unfiltered spectral response of the ROTSE-I 
imaging system which spans a large part of the Johnson- 
Cousins system from mid B to mid / (Section 2.2). The 
best empirical prediction of a ROTSE magnitude TUy^^^g^ 
for Tycho stars is: 

TOg -niy 

fn — nn — — — — 
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The median shift between instrumental magnitudes and 
the above color-corrected magnitudes of Tycho stars is 
then applied to all stellar magnitudes from a given im- 
age. This procedure puts ROTSE-I measurements onto 
a V^-equivalent scale, in the sense that the mean Tycho 
star has Wv-rotse ~ "^v These intermediate object lists 
are passed to light curve building software and are subject 
to further refinements (Section 3.2). Matching to Tycho 
stars occasionally fails in very crowded areas of the Galac- 
tic plane and/or due to substantial cloud cover. Initial 
calibration was successful for 184,006 frames. 

3.2. Object Identification and Photometric Corrections 

A set of object lists for all exposures of a given field has 
to be collated in order to identify measurements that be- 
long to the same objects and hence construct light curves. 
In the process, a collective look at the temporal behav- 
ior of stars across the field can provide useful information 
on systematics of the photometry. Very wide field images 
show some complications that axe usually unimportant in 



data from narrow field instruments. Pronounced gradi- 
ents of background, color dependent atmospheric; extinc- 
tion and gray extinction from thin clouds are common. 
Color dependent effects cannot be fully corrected using 
only a single photometric band. Some ROTSE-I images 
have additional complications near frame edges caused by 
shutter problems (Section 2.2). After the processing steps 
outlined in Section 3.1, most photometric residuals are 
still correlated over spatial scales of a few hundred pixels. 
These gradients are handled by local photometric correc- 
tions. The procedure allows measurements of an object 
originating from each frame to be expressed on a relative 
scale with respect to stars in the neighborhood, as given 
by a master list of stars called a template. At the same 
time we can collect diagnostics providing useful measures 
of the data quality in the final database. They are used 
to set the measurement quality flags (Table 2), making it 
easier to select measurements satisfying requirements of a 
particular application. 

3.2.1. Template Construction 

The process begins with preparation of a template ob- 
ject list for each of the 644 fields. Frames with the 
MOUNTFLIP flag (see Table 2) or fewer than 1000 detected 
objects are rejected for this purpose. We also apply a cut 
on standard deviation of the position and magnitude off- 
sets around the fit to Tycho stars (pos_sigma < 0.3 pixels 
and zp_sigma < 0.4 mag, see Section 4.2). Measm-ements 
of all field objects from all admitted images are clustered in 
order to identify persistent objects. To accomplish this, we 
count the number of individual detections on a simple grid 
with pixels equal in size to detector pixels. Local maxima 
of this histogram containing at least 15% of all possible 
detections are declared as objects. From this counting ex- 
ercise we derive crude centroids that are further refined 
using a 1-pixel identification radius and subsequent sigma 
clipping. The main properties of objects stored in a tem- 
plate list for each field are median (a2000) <^20oo) positions, 
median object magnitudes and standard deviation of in- 
dividual magnitudes around the median. These are essen- 
tially aggregate parameters of the preliminary light curve 
for each object. By using median object magnitudes to 
construct the template, we avoid biasing all field photom- 
etry toward a single frame that may be affected by sys- 
tematic effects of the kind described above. As a measure 
of scatter that is robust against occasional strong outliers, 
we adopt half size of a centered range that includes 68% 
of all magnitude points. 

3.2.2. Photometric Correction Maps 

In the next step, template lists are used to derive 2D 
maps of photometric corrections. The most significant 
contribution to systematics in frame-to-frame photometry 
is caused by clouds, atmospheric extinction gradients and 
occasional shutter problems. The photometric errors are 
therefore correlated over the scales of a few hundred pixels. 
By examining a median shift between magnitudes of stars 
detected in a program frame and template magnitudes for 
the same stars in any 200 x 200 pixel area, one obtains a 
correction to the photometry of stars in the defined macro- 
pixel. Stars for which the scatter is more than 4 times the 
median error bar are rejected as potentially variable. The 
frames are then divided into macro-pixels, each covering 
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200 X 200 detector pixels, to obtain 10 x 10 maps. The map 
records relative photometry corrections, the scatter of all 
magnitude differences in a macro-pixel and the number of 
objects available for those calculations. The global scatter 
of the map is also recorded. All these characteristics can 
be used to assess the quality of the data. The pipeline will 
assign a null value to a macro-pixel in cases when there 
are fewer than 10 stars to work with, available stars cover 
less than half of the macro-pixel area, or the value of the 
correction or its error are unreasonably high (> 1.0 mag). 

3.2.3. Compilation of Light Curve Database and Object 

Catalog 

Applying relative photometry corrections and building 
intermediate database files is the final stage in the data 
reduction process. Individual corrections are obtained by 
bi-linearly interpolating the maps. In the case of a miss- 
ing macro-pixel surrounded by valid ones, we allow a patch 
based on linear interpolation to be applied. Such measure- 
ments are flagged accordingly but were observed to remain 
reliable. We also flag measurements with large corrections 
(> 0.1 mag), large scatter of magnitude differences used 
to derive the correction (> 0.2 mag), and large scatter of 
all macro-pixels in the map (> 0.1 mag). If the patch 
cannot be computed, there is insufficient information to 
derive the correction. This condition sets the NOCORR flag 
signaling a very unreliable measurement. Table 2 describes 
all the data processing flags. The correction coifld not be 
calculated for fewer than 0.2% of the database measure- 
ments. For most frames 95% of all correction values are 
distributed within ±0.05 mag. In case of a frame affected 
by a shutter glitch the distribution can be strongly non- 
Gaussian with the 95% interval typically extending to ±0.1 
mag. In both cases there is a small tail of high correction 
values extending to about ±1.0 mag. 

Measurements that can be associated with the position 
of a template object to within 4a are tagged with the ob- 
ject ID and frame ID and stored in the main light curve 
database file. The error estimate is available from the 
template list. There are no additional restrictions on de- 
tections admitted to the main light curve database. The 
final number of database light curves is almost 2x10'' 
(the same as the number of template objects). There are 
^3.35x10^ individual observations with identified parent 
objects. Remaining observations are in the "orphan" cat- 
egory. We decided to keep only high S /N orphan measure- 
ments, those with magnitudes 14.5 or brighter and magni- 
tude errors below 0.1. The database contains over 2x10® 
orphan measurements. All measurements are treated indi- 
vidually, in particular joined detection in paired exposures 
(Section 2.3) is not required. Such criteria, based solely 
on frame epochs, can always be applied to the data as a 
post processing step. 

At this point we formed a revised catalog of objects with 
average properties based on the fully processed light curve 
files. For each template object we examine a corrected set 
of measurements and recalculate the median magnitude, 
magnitude scatter (significantly improved in most cases), 
median error bar and other useful characteristics. Only 
the subset of "good" measurements defined in Table 3 is 
used in these calculations. These criteria primarily reject 
error codes and several flags associated with known prob- 
lems. Currently, it is a recommended choice for working 



with NSVS data that offers a sensible compromise between 
the amount of data and data quality. Other equally good 
or better selection cuts may be possible. A small fraction 
of all objects below the equator and at the far field edges 
have fewer than 15 "good" measurements and, as such, 
their data are of limited use. Corresponding catalog en- 
tries have light curve statistics copied from the template 
list and are flagged as having TPLSTATS. 

3.2.4. Multiple Detections in Overlap Regions 

Pointing imperfections in repeated exposures of a given 
field result in object detections outside the area around 
the field center covered by a single CCD format (Table 1). 
This effectively extends each template by a margin of ±150 
pixels on all sides. Together with similarly sized intended 
overlap between the fields, this causes multiple detections 
of the same physical objects in more than one field. Over- 
lap regions tend to grow closer to the celestial north pole 
for simple geometric reasons. About 30% of the database 
light curves belong to object references of second order 
or higher. In Section 4.2 we use those independent light 
curves to assess systematic errors in positions and magni- 
tudes. On the other hand, removing multiple references 
to the same objects requires perfect control over system- 
atic effects. Even if the influence of shutter problems and 
atmospheric gradients were known with certainty, some 
irreducible effects would remain. The problem is not triv- 
ial in the case of unfiltered photometry, especially for fast 
variable sources, but it becomes much easier to handle for 
a few sources of particular interest after the data has been 
extracted from the database. This is why a global merge 
of the physical database has not been performed. We pro- 
vide a database table with pairs of object IDs that refer 
to the same physical object (Section 4.2). Implementation 
details of this kind can be hidden by the user interface to 
the database by defining a mapping between the low level 
database structures and the top level view where object 
data appear as merged (Section 4.2). That approach pro- 
vides full flexibility to refine the definition of what is the 
same, and what is not, as more is known about the data 
set. 

4. PUBLIC DATA DISTRIBUTION 

4.1. Survey Coverage and Quality 

Our Northern Sky Variability Survey covers the entire 
sky visible from Los Alamos, New Mexico. This includes 
the entire region north of the declination S ~ —38°; more 
than 30,000 square degrees, 75% of the celestial sphere. 
The completeness, quality of the data and the number of 
available measurements are noticeably lower at low decli- 
nations and low Galactic latitude. The overall quality of 
the survey measurements is summarized in Table 4. In 
Figure 2 we show positions of roughly 5x10^ NSVS stars 
brighter than 1 1 mag out of about 14 million stars brighter 
than 15.5 mag present in the survey. The main feature in 
this plot is the overdensity of stars near the plane of the 
Milky Way. In that region one can notice relatively less 
populated areas due to dust lanes. An empty spot near 
the Galactic Center (field 156 d) is an artifact. The re- 
gion is basically lost from the survey due to the very low 
number of available frames. Each dot in Figure 2 repre- 
sents a time series of typically a few hundred points spread 
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over a total time baseline of 1 year. Temporal coverage is 
subject to yearly visibility patterns. The trade-ofF between 
this enormous monitoring coverage on one side, and on the 
other, relatively low resolution compounded by complex- 
ities of very wide field photometry, shaped the final data 
quality. The faintest objects recorded in the survey have 
V ~15.5 mag, however the incompleteness starts increas- 
ing sharply at 15th magnitude. Saturation may occur in 
stars as faint as 10.5 mag on some nights, but mostly af- 
fects stars brighter than 10 mag. Due to exposure times 
shortened by a factor of 4 during bright time, some stars 
as bright as 8.0 mag still have a number of unsaturated 
measurements sufficient for analysis. 

Figure 3 summarizes important statistics of the NSVS. 
Numerous survey parameters strongly correlate with the 
Galactic latitude and declination. Several artifacts due 
to the low volume and quality of available data near the 
Galactic plane at low declinations are visible. The bright 
spot in the plots of number of available frames is due to 
a special data run in field 072 (Kehoe et al. 2002). The 
darker areas near the Galactic plane are a result of a lower 
success rate in matching frames to the Tycho catalog. 
Strong differences near the celestial pole between panels 
c) and d) are caused by exclusion of flag MOUNTFLIP in the 
definition of a "good" measurement. Lower than average 
performance of camera d is evident from a periodic pattern 
in the map of photometric scatter (panel b). The field pat- 
tern in the number of database frames reveals that camera 
c was not collecting data for about 3 months, when it had 
to be serviced after an electronics failure. Although the 
number of useful measurements below the equator drops 
dramatically to fewer than 100, this is still sufficient to 
detect variables and the NSVS remains useful over a large 
section of the southern hemisphere. 



4.1.1. Astrometric Errors 

Statistical scatter of object positions in individual 
frames can be better than 0.1 pixel (1.4") within a sin- 
gle field (lower part of Figure 4) . Median positions from 
a large number of measurements should be much more ac- 
curate than that, however that ignores systematic errors 
of the coordinate system derived separately in each field. 
A better measure for overall positional accuracy is the dif- 
ference of median positions (in 2-D this time) for the same 
bright unsaturated stars observed in overlapping parts of 
adjacent fields. Such differences should be dominated by 
a systematic contribution and are shown in Figure 5. The 
distribution of these offsets turns out to be comparable to 
that from Figure 4 and fits well within a single image pixel. 
Figure; 5 also shows how typical position uncertainties are 
affected by the high density of stars in the vicinity of the 
Galactic plane. The distribution peaks at a higher value 
for the error and develops a much longer tail due to strong 
blending. 

4.1.2. Photometric Errors 

Figure 4 (upper panel) presents magnitude scatter as 
a function of median object magnitude in a random field 
away from the Galactic equator. Photometric scatter is es- 
timated using "good" photometric points (Section 3.2.3) 



corresponding to about 75% of the best data. We consider 
magnitude scatter for a median star between 11 and 12 
mag in each field to be the "limiting scatter" , that is the 
best attainable for a significant fraction of bright stars in 
a given field. Observations in this magnitude range are 
limited by various systematics of photometric conditions 
and data reductions rather than the statistics of the back- 
ground. There is a strong correlation between limiting 
scatter and Galactic latitude which is evident in Figure 6. 
The median limiting scatter over the total survey area 
is 0.02 mag. However the Galactic plane region, where 
photometric accuracy suffers significant degradation, has 
a higher number density of objects. Averaged over an 
8° X 8° field, spatial density of the NSVS objects can vary 
between 150 and 1000 deg~^. 

Similar to astrometry, systematics of photometry can be 
investigated by examining the offsets between magnitudes 
measured for stars in overlapping parts of adjacent fields. 
Histograms of these differences for two random overlap re- 
gions at low and high Galactic latitudes are dominated 
by constant stars, and are shown in Figure 5. Half of 
the stars in the figure differ by 0.04 mag or less in their 
median magnitude obtained from light curves constructed 
independently in neighboring fields. In individual cases, 
however, suc;h differences can reach 0.1-0.2 mag. They 
result primarily from residual shutter problems that prop- 
agated to the construction of field templates, but also from 
irreducible atmospheric color effects in single broad band 
photometry. The fact that the histogram width does not 
change in a significant way in the proximity of the Galac- 
tic equator agrees with our assessment that the differences 
arise due to systematics associated with a particular set of 
field images and not because of high number density of 
stars or blending. It must be stressed that the overlap re- 
gions between the fields provide the upper bound on the 
systematic errors. This is where many instrumental effects 
(Section 2) manifest themselves strongest. The internal 
consistcincy over the remaining area is certainly much bet- 
ter, although not easily studied without a suitable external 
comparison catalog. 

4.2. Public Database 

All photometric time series data in the NSVS is available 
for public access. The primary means to search and extract 
the data is Sky Database for Objects in Time-Domain 
(SkyDOT^; Wozniak et al. 2002). SkyDOT is intended to 
become a virtual observatory for variable objects. The site 
will provide a uniform interface to several large temporal 
data sets with a number of time series analysis tools for on- 
the-spot application to currently viewed data. SkyDOT 
is implemented using PostgreSQL^, an Object Relational 
Database Management System with support for practi- 
cally the entire SQL92 standard. The present database 
includes six entities listed in Table 5. The columns are 
explained in Tables 6 and 7. Five of those tables represent 
major entities in temporal work: Field, Frame, Object, 
Observation, Orphan. The remaining one, Synonym, helps 
to identify multiple references to the same physical objects 
(Section 3.2.4). It implements "is the same" relationship 
between entries of the Object table. Users familiar with 
SQL can submit queries to the engine with some limita- 
tions on the size of output imposed. All users can access 
the database through a graphical interface offered by the 
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website. The most popular SQL queries can be aceessed 
through browser buttons and search forms. Currently, 
only NSVS photometry is available in the public domain. 
There is neither a browsing capability nor a direct down- 
load option for ^2 TB of the survey image data. However, 
image access is planned for future versions of SkyDOT. 

We take advantage of limited precision and range of 
some database quantities (like magnitudes) and store them 
as short integers rather than floating point numbers. This 
simple rescaling trick generates close to 50% storage sav- 
ings due to a small record size in the main table that dom- 
inates the size of the database. Only users of low and in- 
termediate level data products need to be concerned with 
these technical details, since they are transparent to top 
level users who access data through the presentation layer. 

The data is searchable by celestial coordinates. The Hi- 
erarchical Triangle Mesh (HTM; Kunszt, Szalay & Thakar 
2001) is used for very rapid indexing of positions on the 
sphere. The NSVS database utilizes HTM partitioning 
of depth 14 with corresponding ^1' cell size. Searching 
for position matches within a small tolerance radius can 
be performed at a rate of ~20,000 matches per minute. 
Access time is nearly uniform across the sphere. In partic- 
ular, extraction of objects in regions including the north 
celestial pole is handled seamlessly. HTM runs as a shared 
library extension to the database server with SQL wrap- 
pers. 

4.3. Other Options for Data Access 

All data discussed in the present release can also be 
downloaded from the SkyDOT website'', however transfer 
of the entire database over the network may be impracti- 
cal. Therefore, we will also, upon request, distribute the 
data on user provided DDS 4 tapes (DAT size). The most 
efficient way of storing the database as files in terms of 
the size and speed of rebuilding the database proved to 
be gzip compressed, semicolon separated ASCII files. To 
save space, we skip columns with row IDs running from 1 
to number of entries. Those can always be deduced. The 
total size of the dataset in this form is around 35 GB. 
It is strongly dominated by the Observation table with 
all object light curves, followed by the Orphan table with 
observations that could not be identified with any of the 
objects. Each of the tables Field, Frame and Object fit 
within a single file. Tables Observation and Orphan are 
broken up into 644 files, one for each field. 



4.4. Acknowledging NSVS Data 

We request that the researchers using NSVS in their 
published work include the following statement to ac- 
knowledge the source of data: 

This publication makes use of the data from the North- 
ern Sky Variability Survey created jointly by the Los 
Alamos National Laboratory and University of Michigan. 
The NSVS was funded by the US Department of Energy, 
the National Aeronautics and Space Administration and 
the National Science Foundation. 



5. DISCUSSION AND FUTURE WORK 

We presented the Northern Sky Variability Survey, to 
this date the most extensive temporal record of the sky 
on large spatial scales. All of the survey data is available 
to the astronomical community and can be searched effi- 
ciently using the public SkyDOT database. The database 
contains a total of 3.35 billion measurements for approxi- 
mately 14 million objects in the 8 15.5 magnitude range. 
Time sampling over one full year is between twice per 
night and once every four nights on average. ROTSE-I in- 
strument has achieved a complete spatial coverage of the 
northern hemisphere and a large fraction of the southern 
sky using remarkably low cost hardware. These two fac- 
tors pose limits to the level of detail at which variability of 
the sky was recorded: low spatial resolution, spatial sen- 
sitivity variations, a non-standard filter, and complicated 
systematics near the Galactic plane. 

Despite its limitations, the NSVS is a truly rich source 
of information on stellar variability. Among stars in 
the Galaxy, the fraction of variables with amplitudes de- 
tectable by the NSVS is about -1% (Eyer 1999, Eyer & 
Cuipers 2000). Based on that and on preliminary results 
in Akerlof et al. (2000a), one can expect that tens of 
thousands of new variable stars with good uniform quality 
light curves are present in the data set. Current database 
schema needs to be expanded along the lines described in 
Wozniak et al. (2002) to accommodate various types of 
variables and provide classification capability. The NSVS 
combined with astrometric catalogs providing distances 
and motions, as well as multicolor surveys (2MASS, or 
even SDSS in a narrow magnitude range) will enable a 
comprehensive look at the Galaxy as traced by variable 
stars. The NSVS objects are bright and therefore the pre- 
ferred targets for detailed spectroscopic and astrometric 
work. Spatial resolution of the survey is not far from that 
of high energy sky catalogs like the ROSAT All Sky Sur- 
vey (RASS; Voges et al. 1999) or the XMM Catalog of 
Serendipitous Sources (Watson 2003), and therefore well 
suited for cross correlations. Perhaps the most exciting 
questions to be attacked using the NSVS are about rare, 
hard to find objects. The astronomical literature provides 
numerous unexplained reports of variability events on nor- 
mal stars (e.g. Schaefer, King & Deliyannis 2000 and ref- 
erences therein). 

Photometric monitoring data for Active Galactic Nuclei 
(AGN) providing diagnostics of accretion flows is valuable, 
but limited. Only a few bright AGN are within the mag- 
nitude limit of the NSVS so the real contribution to AGN 
physics will require deeper flux limits and better resolution 
in future projects. A major but low cost improvement in 
data usability would be the use of a set of standard filters 
before starting deeper surveys with more frequent time 
sampling. 

Small robotic telescopes with automated data process- 
ing pipelines are the best candidates for closing the gap in 
the current level of temporal monitoring of the sky. The 
computing power to perform on-line photometry is avail- 
able. Experiments like RAPTOR (Vestrand et al. 2002) 
are starting to tackle the problem of real time detection 



'^http: / / skydot.lanl.gov 
®http; / /www.postgresql.org 
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and immediate follow-up of short time-scale phenomena. 
One can envision a monitoring system capable of par- 
tial interpretation of various events occurring on a variety 
of time-scales and notifying subscribers about interesting 
changes of objects in their scientific problem domain. The 
main challenge is making the immense data stream com- 
prehensible by putting enough smarts into the software. 
The sky itself is the ultimate astronomical database that 
should be mined continuously and in real time. 



This work was supported by the Laboratory Directed 
Research and Development funds at LANL under DOE 
contract W-7405-ENG-36 to the RAPTOR project and 
NASA grant NAG5-5281 to the ROTSE-I collaboration. 
K.K. was supported by the NSF grant AST-0205813 to 
Michigan State University, and S.M. was supported un- 
der the auspices of the DOE, NSSA by UC, LLNL under 
contract W7405-ENG-48. 
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Table 1 

Equipment and Operations in the Northern Sky Variability Survey 



Parameter Value 



Telescope and Site 

Geographic position Los Alamos, New Mexico: W106°15'13", N35°52'9' 

Elevation 2300 m 

Mount speed 100°s-i 

Telescopes Four 200mm, f/1.8 Canon lenses 

Seeing FWHM-20", instrumental 

Vignetting Up to 40% near frame corners 

Lens offsets Aacos((5), A6 with respect to mount position 

a: -4°,+4° 

b: -4°, -4° 

c: +4°, -4° 

d: +4°, +4° 



Imaging Cameras: Four Apogee AP-lOs 



CCDs 2kx 2k Thomson TH7899M chip 

Gain 8 e- ADU"! 

Read noise 13-25 e~ pixel"^ 

Dynamic range > 74 dB 

Read mode 14 bit, 1.3 MHz 

Frame size 2035x2069 pixel 

Image scale 3.50 mm/° 

Pixel size and scale 14(U; 14.4"pixel~^ 

Filter Unfiltered optical response ~450nm-1000nm, 

effective wavelength of R band 

Field of view 8.2° x 8.2° per camera 



Data collection 

Survey area 33,326 deg^ {6 > —38°), best coverage for S > 0° 

Time baseline 1 year, from April 1, 1999 to March 30, 2000 

Number of fields 644=161x4 (cameras abed) 

Number of frames 184,006 

Number of nights 275 out of 365 

Time sampling Pairs of frames 1.5 minutes apart, up to 2 pairs per night 

Calibration 500-1000 Tycho stars per frame 
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Table 2 
Processing Flags 



Hexadecimal Decimal Flag Description 
Bit Bit 



Frame flags 

0x0001 1 MOUNTFLIP Mount flip near north pole, fields ab observed by cameras cd 

0x0002 2 ODDMNTPOS Non-standard mount position, only fields 001 abed in year 2000 



Object flags 

0x0001 1 TPLSTATS Light curve statistics in Object table without photometric corrections 

0x0002 2 BIGSHIFT Final median object centroid more than la from preliminary position 



Measurement Flags: SExtractor 

0x0001 1 NEIGHBORS > 10% of object area affected by a neighboring object or bad pixels 

0x0002 2 BLENDED Object is a result of deblending procedure 

0x0004 4 SATURATED Object has at least 1 saturated pixel 

0x0008 8 ATEDGE Object is tnuicatcd by image boundary 

0x0010 16 APINCOMPL Aperture data incomplete or corrupted 

0x0020 32 ISINCOMPL Isophotal data incomplete or corrupted 

0x0040 64 DBMEMOVR Memory overflow occurred during deblending 

0x0080 128 EXMEMOVR Memory overflow occurred during extraction 



Measurement Flags; Photometric Correction 

0x0100 256 NDCDRR Relative photometry correction could not be calculated 

0x0200 512 PATCH Map of relative photometry corrections was patched to derive correction 

0x0400 1024 LONPTS Low number of points in a macro-pixel (<10) 

0x0800 2048 HISCAT High scatter of magnitude differences in macro-pixel (>0.2 mag) 

0x1000 4096 HICDRR High value of correction (>0.1 mag) 

0x2000 8192 HISIGCORR High scatter of corrections across the map (>0.1 mag) 

0x4000 16384 RADECFLIP Mount flip near north pole occurred (fields 001-032 abed only) 



Table 3 

Definition of a Good Photometric Point 



Condition Serves primarily to remove 



5.0 < magnitude < 16.0 .... SExtractor error codes 

0.0 < mag error < 0.4 SExtractor error codes 

! SATURATED Saturated measurements; important for stars in and out of saturation 

! NOCORR Generally unreliable macro-pixels 

! LONPTS Measurements from far frame edges or shallow images 

! HISCAT Macro-pixels with unreliable photometry 

! HICORR Macro-pixels substantially affected by cloud cover or shutter problem 

! HISIGCORR Frames substantially affected by cloud cover or shutter problem 

! RADECFLIP Possible camera dependent systematics; harmless in many applications 
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Table 4 
NSVS Data Quality 



Parameter Value 



Pointing offsets <t ~ 0.3° (75 pixels) per coordinate 

PSF FWHM^l.S pixels, undcrsamplcd, spatially variable, temporally stable 

Saturation 10-10.5 mag, up to 8 mag in bright time 

Limiting magnitude ~15.5 mag, 14.5 mag in bright time 

Astrometric errors 

Random^ magnitude dependent, 1(7=0.7-4.3", la=1.4-5.8" at Galactic |6| < 20° 

Systematic^ median deviation ~1.2" in general field, ~2.5" at Galactic |6| < 20° 

Photometric errors 

Random'' limiting scatter lcr=0.02 mag in median field, lcr=0.02-0.05 mag at Galactic |6| < 20° 

Systematic'' median deviation ^0.04 mag near frame edges, up to 0.2 mag in extreme cases, 

improving near field center 

Blending Stars closer than 3 pixels (84") generally merged, 

severe at Galactic |6| < 20°, ^5% loss of survey area 

Time sampling Twice per night to once every ~4 nights 

Number of epochs ~200 for average light curve, follows yearly visibility across the sky 

Number of objects ~14 million 



^frame to frame within the same field 

^systematic difference for the same object in overlap region between fields 



Table 5 
Database tables 



Table Number of rows Description 



Field 644 ROTSE-I patrol tiles, each camera counted separately 

Frame 184,006 Image header and post-processing frame quality information 

Object 19,995,106 One record of aggregate information for each light curve in the database 

counting separately same object detections from different fields 

Synonym 14,582,566 Pairs of light curve IDs referring to the same physical object 

Observation 3,353,171,900 All measurements for all light curves 

Orphan 208,106,474 Measurements unidentified with any of the light curves 

but brighter than 14.5 mag and with errors < 0.1 mag 
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Table 6 

Explanation of table columns: tables Field and Frame 



Column Name Data Type Unit Description 



Field Tal)l( 



id 


int32 




Field ID (primary key) 


ncime 


char [4] 




Field name 


rac 


float32 


deg 


Q!20oo of field center 


decc 


float32 


deg 


(^2000 of field center 


glc 


f loat32 


deg 


Galactic I of field center 


gbc 


f loat32 


deg 


Galactic h of field center 


nobs 


int32 




Number of frames 


nobj 


int32 




Number of catalog objects 


sig-ph 


f loat32 


map; 


Liinitinp; photometric scatter 


Frame Table 


id 


int32 




Frame ID (primary key) 


field id 


int32 




Field ID (foreign key) 


fname 


char [20] 




Image file name 


camera 


char [2] 




Camera ID (abed) 


mjd 


f loat64 


day 


JD - 2450000.5 


obstime 


f loat64 


s 


Time of observation in UT seconds 


date obs 


char [22] 




UT date of observation 


exptime 


float32 


s 


Exposure time 


bkg 


f loat32 


. . . counts 


Sky background 


bkg_sigma . . . 


f loat32 


. . . counts 


Standard deviation of sky background 


pos_sigma . . . 


f loat32 


pixels 


Standard deviation around the fit to positions of Tycho stars 


zp_offset ... 


float32 


mag 


Median magnitude offset with respect to Tycho stars 


zp Sigma 


f loat32 


mag 


Standard deviation of magnitude offsets 


m lim 


f loat32 


mag 


Limiting magnitude 


sat_mag 


f loat32 


mag 


Saturation magnitude 


nobj det 


int32 




Number of objects detected by SExtractor 


nob j _ext 


int32 




Number of objects actually measured 


nmatch 


int32 




Number of Tycho stars used in magnitude matching 


dmoon 


f loat32 


deg 


Angular distance between frame center and the Moon 


elev 


float32 


deg 


Elevation of frame center with respect to horizon 


azimuth 


float32 


deg 


Azimuth of frame center 


mount jra 


float32 


deg 


0^2000 position of telescope mount 


mount _dec . . . 


f loat32 


deg 


(52000 position of telescope mount 


offst_ra 


f loat32 


deg 


Aa cos((5) offset between frame center and mount position 


offst_dec ... 


float32 


deg 


AS offset between frame center and mount position 


mapjrms 


f loat32 


mag 


Standard deviation of the photometric correction map 


map_ripix 


intl6 




Number of valid pixels in photometric correction map 


flags 


intl6 




Frame flags 
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Table 7 

Explanation of table columns: tables Object. SYNONY^L Observation and Orphan 



Column Name Data Type Unit Description 



Object Table 



id 


int32 




Object ID (primary key) 


rao 


f loat64 


deg 


Median a2ooo of object centroid 


sigjrao 


f loat64 


deg 


Standard deviation of individual 012000 positions 


deco 


f loat64 


deg 


Median 62000 of object centroid 


sig deco 


float64 


deg 


Standard deviation of individual (520oo positions 


htm id 


int64 




HTM ID for quick spatial queries on the sphere 


mag 


float32 


mag 


Median object magnitude from "good" points 


rmsjnag 


f loat32 


mag 


Standard deviation of "good" points around median 


med err 


float32 


mag 


Median error bar of "good" points 


n obs 


intl6 




Number of "good" points 


n_noflip 


intl6 




Number of all points without RADECFLIP flag 


n points 


intl6 




Number of all object detections 


flags 


intl6 




Object flags 


Synonym Table 


idl 


int32 




First object ID (composite primary key) 


id2 


int32 




Second object ID (composite primary key) 


separation . . 


float32 


deg 


Spherical distance (<1 pixel) 


Observation Table 


obj id 


int32 




Object ID (composite primary key & foreign key) 


frsune id 


int32 




Frame ID (composite primary key & foreign key) 


dra 


intl6 


....1/32000 deg 


02000 offset from median centroid 


ddec 


intl6 


....1/32000 deg 


^2000 offset from median centroid 


mag 


intl6 


mmag 


Corrected object magnitude 


err 


intl6 


mmag 


Magnitude error 


flags 


intl6 




Measurement flags 


Orphan Table 


frame id 


int32 




Frame ID (composite primary key & foreign key) 


rao 


float64 


deg 


O2000 of object (composite primary key) 


deco 


f loat64 


deg 


<^20oo of object (composite primary key) 


mag 


intl6 


mmag 


Corrected object magnitude 


err 


intl6 


mmag 


Magnitude error 


flags 


intl6 




Measurement flags 


htm id 


int64 




HTM ID for quick spatial queries on the sphere 
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raw images, 
calibration frames 
(dari<s and fiats) 



templates for aii fields: median 
object magnitude, position, 
magnitude scatter ^ 



corrjmage 

dark subtraction, flat 
field correction, initial 
position and magnitude 
calibration (Tycho stars) 




mktpl 

cluster detections, 
find density peal<s, 
apply Sigma clipping, 
get one entry per object 




apply photom. corrections, set quality flags, 
make final identification of measurements to 
objects, liandle unidentified detections 



object catalog: one entry 
of aggregate info per object 



relphot 

find mag offsets core- 
lated over 200x200 pix 
regions between frames 
and template 



mkcat 

derive aggregate info 
for eacti object using 
good pliotometry, 
set additional flags 



10x10 maps of photometric 
corrections for all frames 



Fig. 1. — The NSVS data processing pipeline. 



This figure is available in gif format: figure2.gif 



Fig. 2. — Positions of NSVS objects brighter than 11 mag (about 500,000) in equal area MoUweide's projection. Each dot represents an 
NSVS object with a temporal history typically composed of a few hundred measurements covering the 1 year baseline. 
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This figure is available in gif format: figure3.gif 



Fig. 3. — NSVS at a glance. Four panels show all sky gray scale maps smoothed over 8° spatial scale; a) Number density of NSVS objects, 
b) photometric scatter for bright unsaturated stars calculated using "good" measurements, c) median number of points per light curve, d) 
median number of "good" photometric points per light curve. 



This figure is available in gif format: figure4.gif 



Fig. 4. — Random errors as given by frame to frame scatter in a typical field. Photometric errors (upper) and position errors (lower) are 
shown as a function of median object magnitude. 
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position difference (deg) mognitude difference 




position difference (deg) mognitude difference 



Fig. 5. — Systematic errors as given by differences between multiple detections of the same objects in overlap regions between the adjacent 
fields in the Galactic plane (lower) and near the Galactic pole (upper). Shown are differences in median object positions (left) and median 
object magnitudes (right) for bright unsaturated stars. Position differences fit well within a small fraction of a pixel. Magnitude offsets result 
from residual shutter problems and properties of very broad band photometry in the presence of intrinsic spread of object colors. Estimates 
based on overlap regions, where numerous instrumental effects are strongest, provide an upper bound on systematic errors. 
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Fig. 6. — Limiting photometric scatter as a function of the Galactic latitude. In each field we show light curve scatter for bright unsaturated 
stars calculated using "good" measurements. 
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